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SUMMARY 

A  description  is  given  of  the  process  used  to  develop  and  verify  the  aerodynamic  data¬ 
base  of  the  F-lllC  flight  dynamic  computer  model.  The  process  uses  stability  and  control 
derivatives  obtained  from  flight  trials  performed  at  the  RAAF's  Aircraft  Research  and 
Development  Unit  (ARDU).  Model  response  has  been  verified  against  aircraft  time  histories 
measured  during  the  trials.  The  results  presented  compare  the  aircraft  response  in  lateral 
manoeuvres,  at  various  wing  sweeps.  A  very  good  degree  of  matching  is  possible,  although 
further  investigation  should  correct  inaccuracies  which  can  cause  poor  estimation  of  yaw 
rate. 
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Introduction 


A  >:ix  degree  of  freedom  flight  dynamic  model  of  the  F  lllC  aircraft  has 
been  developed  at  the  Aeronautical  Research  Laboratory  ( ARL)  to  support 
aircraft  operations  and  aircraft  and  weapon  system  development  (Reference 
1).  The  model  has  been  developed  from  the  manufacturers  design  reports 
and  wind-tunnel  reports  for  the  F-lllA  version.  A  series  of  flight  trials 
was  performed  on  the  the  F-lllC  at  the  RAAF’s  Aircraft  Research  and 
Development  Unit  (ARDU)  to  obtain  stability  and  control  derivatives  for 
validation  of  the  ARL  model. 

Analysis  of  the  data  was  conducted  by  the  Aircraft  Flight  Dynamics 
Group  at  ARL  using  the  parameter  estimation  techniques  described  in  Ref¬ 
erence  2.  Results  from  the  analysis  are  reported  in  Refertnces  3  to  8. 
Small  disturbance  models  of  the  aircraft  longitudinal  and  lateral  aerody¬ 
namic  characteristics  ars  used  in  the  identification  process  to  determine  the 
aerodynamic  parameters.  These  parameters  are  then  used  to  modify  the 
aerodynamic  data  base  in  the  ARL  six  degree  of  freedom  non-linear  flight 
dynamic  model. 

Modifications  have  been  made  to  the  ARL  flight  dynamic  model  to  fa¬ 
cilitate  this  updating  process,  often  referred  to  as  model  validation.  The 
process  is  a  necessary  part  of  flight  dynamic  model  development  br.t  is  not 
rigorously  defined  and  so  relies,  to  a  degree,  on  engineering  judgement.  The 
objective  with  the  F-lllC  model  is  to  combine  the  modified  aerodynamic 
data  base  with  the  existing  model  of  the  flight  control  system  to  give  a  good 
representation  of  the  measured  flight  behaviour. 

The  purpose  of  this  report  is  to  document  this  process  and  to  present 
preliminary  results.  Control  inputs  measured  during  the  flight  tritds  were 
used  to  drive  the  flight  dynamic  model  and  comparisons  were  made  between 
the  time  histories  generated  by  the  model  and  those  measured  in  flight.  This 
report  describes  the  approach  used  for  the  lateral  directional  aerodynamic 
parameters  and  describes  the  evaluation  of  a  number  of  lateral  responses 
from  the  flight  Trials.  A  brief  investigation  of  the  longitudinal  motion  is 
also  included,  which  indicates  that  a  similar  approach  can  be  used  for  the 
longitudinal  parameters. 

2  Model  Description 

The  six  degree  of  freedom  flight  dynamic  model  uses  standard  rigid  body 
aircraft  equations  of  motion.  Quaternions  are  used  to  calculate  aircraft  atti¬ 
tudes  and  a  fourth  order  Runge  Kutta  method  is  used  to  integrate  the  equa¬ 
tions.  The  force  equations  are  calculated  in  airpath  axes  and  the  moment 
equations  are  calculated  in  body  axes.  These  equations  are  programmed 
using  the  Advanced  Continuous  Simulation  Language  (AC(SL)  see  Ref.  9. 
and  are  described  in  detail  in  Ref.  10. 

The  aerodynamic  database  contains  the  stability  and  control  derivatives 
that  allow  calculation  of  forces  and  moments  acting  on  the  aircraft  due  to 
its  own  motion  and  to  control  surface  deflections. 
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3  Procedures  for  Verifying  the  Lateral  Deriva¬ 
tives 

The  lateral  response  of  an  aircraft  may  be  assessed  by  inspection  of  the 
response  variables:  roll  rate  p,  yaw  rate  r.  and  angle  of  sideslip  J.  In 
verifying  the  influence  of  the  lateral  derivatives,  the  calculated  responses 
of  these  variables  to  measured  control  inputs  have  been  compared  with  the 
measured  response. 

The  derivatives  which  influence  the  lateral  directional  response  are  con¬ 
tained  in  the  following  force  and  moment  equations  which  are  used  in  the 
flight  dynamic  model. 

Sideforce  Equation 


FVS  =  Cy^P  +  Cy^r  +  Cy^S  +  Cy,Ja  +  Cy^^tlT  +  Cy,jsp 

^-weight  x  sin  0 cos 0  (1) 

Rolling  Moment  Equation 

i  =  OpP  +  Ct,r  +  C’lgd  +  Ci,Ja  +  C'lfjT  +  C'li^Jsp  (2) 

Yawing  Moment  Equation 

=  Cn,p  +  c„,r  +  Cngli  +  C„,Ja  +  C„,Jt  +  C„„^6sp  |3) 

These  total  forces  and  moments  are  substituted  into  the  flight  dynamic 
model  and  the  state  equations  are  integrated  to  produce  the  output  re¬ 
sponses  in  roll,  yaw  and  sideslip. 

Because  of  coupling  between  all  three  equations  it  is  a  complex  problem 
to  isolate  the  source  of  any  errors.  For  example,  if  the  calculated  rolling 
moment  is  poorly  matched,  tmd  is  traced  to  an  incorrect  value  for  the  roll 
due  to  yaw,  this  could  be  due  to  an  error  in  either  the  value  of  the  derivative 
C'l^  or  due  to  errors  in  the  response  variable  r. 

Some  means  of  isolating  the  calculation  of  the  three  motions  p.  r,  and  0 
was  considered  desirable,  so  that  each  equation  could  be  analysed  separately, 
using  flight  measured  values  for  p,  r  and  0. 

This  has  been  organised  in  the  ACSL  simulation  programme  by  using 
three  logical  switches  called  BEBM  (Beta  Equals  Beta  Measured)  PEPM. 
and  RERM  for  p  and  r  respectively.  By  setting  any  two  of  these  switches  to 
.TRUE.,  and  the  third  to  .FALSE.,  the  equations  are  calculated  using  the 
actual  measured  values  of  the  first  two  variables,  and  the  computed  value 
for  the  third. 

Further  examination  of  the  full  six  degree  of  freedom  equations  of  motion 
reveals  that  there  can  also  be  significant  coupling  from  the  longitudinal  mo¬ 
tion  into  the  lateral  response.  Two  additional  logic  switches  have  therefore 
been  added  for  alpha  and  q,  these  are;  AEAM  and  QEQM. 

Thus  by  setting  RERM=.T.  and  BEBM=.T.,  but  PEPM=.F.,  in  addi¬ 
tion  to  AEAM=.T.  and  QEQM=.T.,  the  roll  rate  p  will  be  calculated  using 
measured  values  of  r,  0,  a,  and  q,  and  so  the  Cj  derivatives  can  be  checked 
without  errors  caused  by  incorrect  calculation  of  these  parameters  affecting 
the  accuracy  of  the  calculated  roll  rate. 
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When  flight  measured  data  is  not  used,  the  Ingiral  flag  FLTDAI  i'  ^<’1 
to  False.  Thus,  having  verified  the  f'„,  f’;.  C\  derivatives  using  selerleil 
measured  responses,  the  niivdel  can  then  he  run  with  all  three  variables 
calculated  simultaneously. 

The  procedure  adopted  to  verify  the  derivatives  is  as  follows: 

1.  Run  the  AC'SL  model  using  database  values. 

2.  Run  model  using  Flight  Trial  values  with  AEAM=.T.  QEQM  =  .T.  and: 

(a|  PEPM=.F.  RERM=.T.  BEBM=.T.,  observe  accuracy  of  p 

(b)  PEPM=.T.  RERM  =  .F.  BEBM=.T.,  observe  accuracy  of  r 

(c)  PEPM=.T.  RERM=.T.  PEBM=.F.,  observe  accuracy  of  J 

(d)  Run  with  PEPM,RERM.BEBM=.F.,  so  p,r,  and  /3  are  all  calcu¬ 
lated  simultaneously. 

Coupling  from  the  longitudinal  motion  into  the  lateral  directional  re¬ 
sponse  can  occur  in  two  ways.  Firstly,  the  lateral  aerodynamic  data  are 
stored  in  the  model  in  derivative  form  as  functions  of  the  Mach  no.,  altitude, 
and  the  longitudinal  variable,  angle  of  attack.  Secondly,  for  manoeuvres  with 
large  rates  of  yaw  and  particularly  roll,  the  sideforce  equation  includes  a  sig¬ 
nificant  inertial  contribution  which  is  a  function  of  angle  of  attack.  The  first 
contribution  is  generally  small  unless  the  longitudinal  variables  are  grossly 
in  error,  however  the  second  contribution  can  be  large. 

For  the  Airpath  system  of  axes  used  in  the  F-lllC  ACSL  progranmie 
the  sideforce  equation  can  be  arranged  as  follows: 

=  FYW  X  rmasslvtas  -  RSB  (4) 

where  FYW  =  FYScos{0)  -  FXS sm(0) 

FXS  =  {xthrust  -t-  weight  sin(S))cos(a) 

+{zthrust  +  weight  cos(S))sin(a)  -  drag 
FYS  =  -f-  Cy^P  +  Cyg0  +  Cy,Ja  A  Cy,jT  +  + 

weight  x  sin^costf 
and  RSB-=  psin(a) -F  r  eos(Q) 


FXS  is  the  force  in  the  X  direction  which  is  calculated  from  the  longi¬ 
tudinal  equations,  and  FYS  is  the  force  in  the  Y  direction  and  is  calculated 
from  the  sideforce  equation  (Equation  1). 

The  RSB  term  has  a  very  significant  influence  on  j  in  high  roll  rate 
manoeuvres,  and  so  will  be  sensitive  to  the  accuracy  of  the  predicted  angle 
of  attack  a. 

The  roll  and  yaw  rates  are  derived  from  the  anguleir  accelerations  p  and 
r,  which  are  calculated  from  the  rolling  and  yawing  moment  equations  as 
follows: 


p=  LCi  +  NC7  +  (K's  +  rO  )g  (5) 

r  =  NCe  +  LC2  +  {pCf  -  rC3)g  (6) 

where  Cx  are  inertial  terras  as  given  by: 


(3) 


(  0 

/.V.V 

hz 

/'vz 

("■l  ^ 

=  Izz 

Co 

('2  = 

=  l.\z 

:c„ 

C.S  - 

-  f  2  ^ 

(/.V.V 

/it 

•  Izz) 

4 

Cl  ■ 

l/v) 

hz) 

(C2 

■  Ixz) 

^  ' 

=  /.V.Y 

,  (  '0 

^'9  ' 

=  f's  / 

(/.V.V 

lyy ) 

*  (C? 

'  /.vzl 

It  can  be  seen  that  g,  the  pitch  rate,  contributes  gyroscopic  inoinents. 
A  switch  QEQM  is  defined  in  the  progranune.  which  when  set  to  .TRI  E, 
substitutes  the  measured  value  of  q  in  these  equations.  Any  differences 
between  the  measured  and  calculated  roll  rate  p  or  yaw  rate  r  can  then 
be  traced  to  either  the  lateral  derivatives  or  to  errors  in  calculation  of  the 
longitudinal  motion. 

4  Implementation  of  Flight  Data  in  Flight  Dy¬ 
namic  Model 

Additional  subroutines  were  written  to  enable  tbe  flight  dynamic  model  to 
access  the  aerodynamic  derivatives  obtained  from  flight  nteasurements. 

For  the  verification  phase  it  is  proposed  to  transform  the  existing  database 
values,  which  are  based  on  wind-tunnel  measurements,  into  the  new  values 
by  applying  a  scale  factor  and  a  bias  terms  according  to  the  equation: 

Cl.^.yNEW  =  Ci,„  yOLD  scale  factor  -!-  bias 

Appropriate  scale  factors  and  bias  errors  were  calculated  as  part  of  the 
flight  test  analysis  nrocedures.  An  alternative  optic-e  nl'o  provided  whir'’ 
switches  the  data  directly  from  the  existing  data  to  the  new  data  by  setting 
the  scale  factor  to  zero  and  the  bias  term  equal  to  the  new  data. 

To  control  this  operation  an  ACSL  logic  parameter  LATFDR  has  been 
defined,  which  if  set  to  .TRUE,  will  then  modify  the  LATeral  Flight  DeRiva- 
tives  to  the  flight  derived  values. 

With  r,.ATFDR  =  .TRUE.  subroutine  LAT.F  will  call  LATDB.F  to  obtain 
the  database  values,  then  will  call  subroutine  MODLATDERV.F  to  MODily 
the  LATeral  DERViatives.  A  selection  of  the  derivatives  to  be  modified  can 
be  made  via  the  subroutine  MODLATDERV.F.  The  routine  reads  two  files. 
MODDERV.DAT  and  LATCOEF.S.  File  LATC’OEFS  contains  the  scale  fac¬ 
tor  and  bias  values  to  be  applied,  whilst  the  MODDERV'.DAT  file  contains 
a  list  of  lateral  derivatives,  and  assigns  a  logic  value  to  each  derviative.  Only 
if  this  logic  value  is  equal  to  .TRUE,  will  the  specific  lateral  derivative  be 
modified. 

Therefore,  to  convert  from  the  existing  data  base  to  the  flight  derived 
data  both  the  ACSL  parameter  LATFDR.  and  also  the  logic  variable  in  the 
MODDERV.DAT  file  must  be  set  to  TRUE. 

Similar  procedures  are  provided  for  the  longitudinal  derivatives,  via  sub¬ 
routines  LONG.F.  LONGDB.F,  MODDERV.DAT  and  LONGCOEFS. 


(41 


5 


Modifications  to  ACSL  Model 


Purina  til**  Vitrification  *‘xrrfisp.  a  niinif>ar  of  r*trr»Tt ion^.  riiariKc-.  .on]  .-ii 
Slant  "lOPiit  werp  niadp  to  the  At  'SI.  iiiu<|pi.  TIip'p  ar*-  '■iiminari'cil  nion  r 
tSip  foliowing  spvpii  (lifferpiil  hpatlings. 

1.  Corrertioii  of  errors.  When  the  Highl  trials  stahility  ami  miiirol 
iSerivalivPs  were  initialv  used,  poor  matching  resulted.  After  investiga 
tion.  a  number  of  errors  were  fouml  in  the  coiling,  which  were  corrected 
a.s  follows. 

I  a)  The  databa.se  contained  values  of  zero  for  the  spoiler  derivatives 
at  certain  spoiler  deflections.  The  values  given  in  the  database 
FU1LATDB.F18  for  '.SPOILER’-4.5  at  ‘SWEEP'- Iv  4r.  were 
zero  for  the  DLSP.  DNSP  and  DV'SP  derivatives.  This  error  was 
rectified  by  replacing  the  zero  values  with  the  values  given  for 
■SPOILER' -26. 

|h)  Negative  spoiler  deflections  were  being  made  positive.  Iti  pro 
gram  LATDB.F  the  spoiler  deflection  value  SPOILER!'  was  tnaile 
an  absolute  value  for  the  pttrpose  of  retrieving  the  vahies  of  the 
derivatives  from  the  database,  but  these  values  were  not  then 
modified  to  the  correct  sign. 

(c)  Limiting  of  the  maxinmm  alpha  value  was  prograimiied  incor 
rectly.  A  minor  error  had  meant  that  the  limiting  would  never 
occur. 

|d)  Centre  of  gravity  corrections  were  being  carried  out  inrorrertly. 
Program  LAT.F  calls  CGLAT  to  apply  eg  corrections,  however 
flight  trial  data  is  corrected  for  instrumentation  offset  in  a  flight 
data  processing  program,  and  so  the  use  of  C'CLAT  was  not  re¬ 
quired  for  the  flight  trial  data.  When  comparing  model  response 
ami  measured  data  it  is  essential  to  make  the  comparison  at  a 
common  reference  point. 

(e|  The  lateral  derivatives  were  being  incorrectly  modified  by  sub 
routine  CONLAT.  This  is  used  for  certain  database  values  in  the 
original  model  development  to  convert  Fill  A  daia  lo  thc..e  of 
an  F-lllC.  but  should  not  be  done  in  the  rase  of  flight  trials 
values. 

(f)  The  variable  break  point  values  for  parameter  MACHI  in  the 
FU1LATDB.F18  database  were  not  being  read  correctly.  These 
values  are  written  on  two  lines,  but  the  second  line  was  not  being 
read. 

2.  Correcting  Control  Inputs.  To  simulate  lateral  manoeuvres  the 
ACSL  program  is  driven  by  the  control  inputs  aileron,  spoiler,  and 
rudder.  In  the  initial  trim  condition  prior  to  the  manoeuvre  the  model 
values  should  aU  be  zero,  since  the  model  trims  in  symmetric  flight. 
However,  in  reality  the  measured  values  will  have  initial  non-zero  val¬ 
ues  as  a  result  of  aircraft  asynunetries  and  biases  in  the  instrumenta¬ 
tion.  This  was  corrected  by  subtracting  the  initial  value  of  each  control 
setting  from  the  control  time  history  records. 
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.5.  Plotting  of  Dorivaf iv«*s.  fh**  arul  a<  r<'<|vn.inin 

'U*rjvativf*s  vv*»r<*  »<»  l>f  Ai  SL  })araiii#‘t*«r'‘.  lhi<  allows  Uihjh  in 

\tv‘  a  V\nn‘tion  of  tnn»'  to  r\i»Mk  tor  m  itw  'l,ii,i 

<'Tr. 

1.  Operation  of  LATFDR,  LONFDR.  riifsf*  A*  SL  lo^ic  paraiiit' 
t»*rs  are  flrtiiiecl  in  the  main  pro>»raiii.  ami  passed  tt>  i  [i*>  1,  A  I  I  aiui 
L()N(i.F  ^nhroutines.  Iftheir  valuer  are  f rue.  then  MODL  AM  >I .H\  I 
or  MODLONCiDERV.F  are  callecl  to  modify  the  lateral  or  longit miinai 
derivatives  to  flight  trial  values. 

■>.  Reading  file  MODDERV.DAT.  This  tile  contains  the  individual 
logic  values  for  the  lateral  and  longitudinal  derivatives.  These  valiu"^ 
are  read  hy  MODLATDERV.F  or  MODLONGDERV.E.  to  ilefennine 
which  of  the  derivatives  are  to  be  inodihed  to  flight  trial  values. 

b.  Reading  file  LATCOEFS,  LONGCOEFS.  Once  program  MOD 
LATDERV.F  has  rear)  file  MODDERV.DAT.  it  then  reads  LATCOEFS 
to  obtain  the  values  of  the  scale  and  bias  factors.  The  use  of  external 
files  such  as  MODDERV.DAT  and  LATCOEFS  allows  these  values  to 
he  changed  without  the  need  to  recompile  and  execute  the  program. 

7.  Operation  of  AEAM.  QEQM.  PEPM.  RERM,  BEBM.  These 
are  logic  parameters  which  have  been  defined  in  the  main  AC.SL  pro 
gram,  and  if  set  to  .TRUE,  will  cause  the  calculated  value  to  be  re 
placed  with  the  tueasured  value. 

An  investigation  was  also  made  into  the  method  used  in  the  model  for 
calculating  spoiler  deflections.  During  the  flight  trials  the  spoiler  deflections 
were  measured  directly.  When  analysing  flight  manoeuvres  these  measure 
ments  are  used  as  inputs  to  the  ACSL  model.  The  ACSL  moiiel  also  cai 
culates  a  spoiler  deflection  from  the  measured  lateral  stick  position.  The 
model  can  thus  use  either  of  these  sources  of  ^.sp  when  simulating  flight  trial 
manoeuvres.  The  best  roll  response  has  been  obtained  using  the  values  of 
>,sp  calculated  from  the  lateral  stick  position,  which  are  in  general  higher 
than  those  which  were  measured  directly.  The  cause  of  this  discrepancy  has 
not  been  identified. 

6  Results 

Comparisons  are  presented  of  the  calculated  response  of  p.  r  and  .1.  using 
flight  trial  lateral  stability  and  control  derisatives.  with  the  measured  re 
spouse.  Figures  1,  2.  .1,  4  and  5  are  shown  for  the  wing  sweep  angles  72.0. 
•'iO.  li.'i,  26  and  16  degrees  respectively.  The  figures  on  the  left  hand  side 
show  p.  r  and  J  when  calculated  separately  using  measured  vahies  of  the 
responses,  whilst  the  figures  on  the  right  show  the  response  for  the  complete 
model.  Prior  to  calculating  the  complete  model  responses,  some  of  the  f  „ 
derivatives  were  altered  to  improve  the  models  estimated  response.  These 
alterations  were  based  on  the  matching  of  yaw  rate  with  the  other  variables 
fixed,  as  shown  (Fig,  6). 

Results  of  a  simulated  longitudinal  manoeuvre  are  shown  iPig.  7).  In 
this  manoeuvre  however  only  database  derivatives  have  been  used.  Again, 
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I  h»'  i^raph^;  on  ilie  Ipf't  hainl  •'how  !h<‘  n'^pon^f*  ii'inu  m.-a'iir*’'!  ..iIih  ' 
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Simulation  of  a  I>utrh  roll  iiianoeuvrf  was  also  pfrlorimul  M  itt.  In 
'tmilar  fashion  to  th»*  longitinlinal  r**siiUs.  onlv  «laiahasf  ihTivativf^  w^*r'• 
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7  Discussion 

7.1  Lateral  Manoeuvres 

l.atfr.il  roll  iiianoPiivres  were  ^inuilateil  for  wing  swpeps  of  Iti.  '2f>.  'lO 
and  72.0  degrras.  Simulalions  were  conducted  vising  ilerivatives  from  the 
databa.se.  and  also  from  flight  trials  data.  Results  are  presenleil  for  ra'^es 
with  i>.  r.  and  .1  obtained  separately  from  flight  or  from  computation. 

For  some  of  these  cases  the  calculated  values  of  p.  r  and  J  are  very  close 
to  the  measured  response  (Fig.  9).  In  these  cases  it  can  be  as.snmed  that 
the  ().  Cn.  Cy  derivatives  are  correct,  and  that  the  ACSL  programs  and 
evpiations  ji  motion  correctly  predicted  the  aircraft  motion.  The  agreement 
also  confirms  that  the  lateral  derivatives  have  been  correctly  estimated  from 
the  flight  data. 

For  some  cases  the  yaw  rate  response  is  incorrectly  calculated  (Fig.  10). 
When  the  three  parameters  p.  r.  and  J.  are  calculated  separately,  using  flight 
trial  derivatives,  the  p  and  J  response  are  very  close  to  measured  behaviour, 
whilst  the  yaw  rate  r,  although  foUowing  the  same  pattern  as  the  measured 
yaw  rate,  does  not  attain  the  same  values. 

The  exact  reason  for  the  poor  accucuracy  of  the  yaw  rate  has  not  been 
established.  Although  altering  the  value  of  the  C'ng  and  f’n,,  derivatives 
can  improve  the  accuracy  of  the  yaw  rate,  when  calculated  in  isolation  (Fig. 
11).  the  matching  deteriorates  again  when  all  three  motions  are  calculated 
simultaneously  (Fig.  12).  The  fact  that  changing  the  C„g  and  deriva¬ 
tive  appears  to  give  better  matching  does  not  prove  these  derivatives  are 
incorrect. 

The  incorrect  calculation  of  the  yaw  rate  also  has  a  degrading  influence 
on  the  accuracy  of  the  other  parameters,  particularly  sideslip  J,  when  all 
three  are  calculated  simultaneously  (Fig.  13). 

The  incorrect  yaw  rate  has  an  effect  on  the  sideslip  due  to  a  coupling 
of  the  motions  as  shown  by  the  equations  used  to  calculate  d.  Whilst 
is  often  a  dominant  derivative,  in  terms  of  it's  numerical  value,  an  investi¬ 
gation  made  by  altering  the  C],  derivative  has  shown  that  ni  Ifii  mpiJ  mil 
manoeuirres  it  has  very  little  influence  on  the  sideslip. 

Investigations  described  in  .Section  3  have  shown  that  in  the  rapid  roll 
manoeuvres  the  RSB  term  has  a  major  influence  on  the  calculation  of  the 
sideslip  angle.  This  is  demonstrated  graphicaly.  by  plotting  the  individual 
terms  which  combine  to  produce  the  sideslip  angle.  As  shown  in  Fig.  14. 
the  calculated  d  (BETAD)  is  very  close  to  the  measured  d  (BETAM).  This 
calculated  d  is  integrated  from  d  (DBETAR)  which  has  two  terms,  one  from 
the  RSB  component,  and  the  other  from  the  FYW  component  (Eqn.  4). 

These  terms  are  plotted  (Fig.  15)  and.  as  ran  be  seen,  the  RSB  term 
is  almost  identical  to  the  negative  of  DBETAR.  and  the  FYW  term  is  very 
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siiiaU  in  comparison  to  J.  This  shows  that,  lor  llu*  lateral  roll  manoeuvre', 
t  is  almost  entirely  caused  by  the  RSB  term  and  not  from  the  FVVV  term 

The  FVW  term  itself  has  two  components,  FVS  and  FXS.  although 
the  F.\S  term  is  very  small  (Fig.  16).  The  FY.S  term  is  calculated  (Fain. 
1 1  from  the  components  due  to  the  derivatives,  and  a  weight  term,  which 
is  dominant.  When  the  C'y  derivatives  are  set  to  zero,  there  is  very  little 
change  in  the  value  of  FYS  calculated  (Fig.  17a),  and  hence  no  detectable 
difference  in  the  J  calculated  (Fig.  17b). 

.Since  the  RSB  term  is  the  dominant  cause  of  the  sideslip  in  the  lateral 
manoeuvres,  then  when  f3  is  calculated  separately,  with  PEPM.  RERM  and 
AEAM  =  .T..  then  the  calculated  J  would  be  expected  to  be  almost  identical 
to  the  measured  J. 

Any  error  in  the  J  angle  calculated  simultaneously  with  p  and  r  is  due 
to  inaccuracies  in  these  values,  and  not  due  to  any  possible  error  of  the  C'y 
derivatives,  which  have  almost  no  effect  for  the  rapid  roll  manoeuvres. 

The  incorrectly  predicted  yaw  rate  could  also  have  an  effect  on  the  roll 
rate  via  the  C'l,  term  used  in  calculating  L  and  the  inertial  term  rC'4  in  the 
equation 


LCi  +  NC2  +  {pC'i  F  rC4)q 

from  which  p  is  calculated.  The  LC’i  term  is  dominant  in  calculating  p. 
since  the  other  components  iVCj  and  (pf's  +  rfTjq  are  effectively  zero.  This 
is  shown  in  (Fig.  18),  where  these  components  are  denoted  if'l,  yC'2  and 
PPRQ  respectively.  This  demonstrates  that  only  the  equation  determining 
L  governs  the  motion  of  the  roll,  and  that  the  only  coupling  arises  from 
the  terms  within  this  equation.  The  accuracy  with  which  the  roll  rate  p  is 
predicted  implies  that  the  rolling  moment  derivatives  are  quite  accurately 
estimated. 

Whilst  the  equation  used  to  calculate  the  yaw  rate  has  a  very  similiar 
form  to  that  used  for  the  roll  rate,  the  yaw  rate  is  often  poorly  predicted. 

The  yaw  rate  acceleration  is  given  by 


r  —  iVCg  +  LC'2  T  {p^ 9  ~ 

and  Fig.  19  shows  the  contributions  to  f  by  each  term.  Unlike  the 
behaviour  of  p,  r  has  contributions  from  all  three  terms  and  has  coupling 
from  the  roll  motion  in  three  terms;  firstly  from  x  p  in  calculating  A’: 
secondly  from  C'l^  x  p  in  calculating  L;  and  finally  the  P  x  f'g  »  term. 

Since  the  roll  derivatives  have  been  shown  to  be  estimated  with  reason¬ 
able  accuracy,  then  the  poor  predition  of  yaw  rate  would  indicate  that  the 
yaw  derivatives  are  not  well  estimated. 

Once  the  yaw  rate  has  caused  inaccuracies  in  the  roll  and  sideslip  values, 
these  incorrect  parameters  will  then  affect  the  yaw  rate,  via  a  feedback  of 
the  errors  in  the  roll  and  sideslip.  This  naturally  complicates  the  system, 
and  makes  it  difficult  to  assertain  cause  and  effect  between  parameters. 

As  shown  in  Fig.  11,  altering  the  C„  derivative  changes  the  yaw  rate 
response,  but  none  of  the  changes  introduced  has  resulted  in  a  satisfactory 
match  for  this  case  (Fig.  12). 

In  some  cases  (Fig.  20)  it  was  noticed  that  the  calculated  response 
shows  a  small  roll  in  the  opposite  direction  to  the  intended  roll  response. 
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1  his  was  shown  to  be  due  to  the  ( term  sinre  all  other  roll  derivatives,  at 
this  point  of  the  manoeuvre,  are  acting  in  the  correct  direction,  whilst  the 
moment  produced  by  »  br  was  the  only  term  opposing  them.  Beiluciiig 
the  value  off'/,,  to  the  database  value  corrected  this  problem  (Fig.  '21). 

The  T'/,,  derivative  may  also  be  in  error  in  other  manoeuvres  hut  is  not 
detectable  since  for  these  cases  the  time  between  initiation  of  rudder  and 
aileron  deflection  is  too  small  to  enable  the  response  to  be  observed.  This 
example  shows  the  complexity  of  the  verification  process.  The  general  policy 
in  these  situations  is  to  accept  the  flight  data,  which  represents  an  average 
of  the  results  over  the  appropriate  Mach  No.  range  for  that  configuration, 

7.2  Longitudinal  Manoeuvres 

An  exploratory  investigation  of  the  longitudinal  model  was  carried  out  using 
the  procedures  developed  in  the  lateral  verification  study.  The  longitudi¬ 
nal  comparisons  shown  in  Fig.  7  were  obtained  using  derivatives  from  the 
database  and  have  not  been  adjusted  to  reflect  flight  trials  results.  However 
it  can  be  seen  that  the  calculated  values  of  9,  o,  and  9  match  up  reasonably 
well  with  the  measured  responses. 

Using  a  similar  procedure  to  the  lateral  analysis,  the  longitudinal  simu¬ 
lations  were  conducted,  with  the  variables  PEPM,  RERM  and  BEBM  set  to 
TRUE.  Although  the  equations  for  calculating  q,  a,  and  9  include  coupling 
from  the  lateral  motion  in  p,  r  and  d,  these  terms  are  generally  small.  How¬ 
ever,  the  facility  has  been  provided  to  enable  measured  lateral  quantities  to 
be  used  if  required. 

Correct  matching  between  the  calculated  and  measured  a  will  depend  on 
the  requirement  that  the  model  trimmed  q.  at  the  start  of  the  simulation,  is 
the  same  as  the  measured  trimmed  a  at  the  start  of  the  manoeuvre.  Part  of 
the  verification  process  will  be  to  establish  a  matching  of  the  trim  conditions. 

7.3  Dutch  Roll  Manoeuvre 

The  estimated  behaviour  of  the  aircraft  in  a  Dutch  roll  is  presented  in  Fig. 
8,  which  shows  reasonable  accuracy  for  p  and  l3  when  they  are  calculated 
in  isolation.  Unfortunately,  when  all  three  variables  are  predicted,  the  fre¬ 
quency  of  the  model  behaviour  is  no  longer  accurate.  These  preliminary 
results  were  obtained  using  F-lllC  database  values.  Better  reults  should 
be  achieved  when  lateral  derivatives  from  the  flight  trials  are  used,  although 
attempts  to  use  these  in  the  Dutch  roll  have  not  been  successful  yet. 
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8  Conclusion 


Procedures  for  updating  the  ARI.  F-lllC  High!  dynaiiiir  iiioilel  vvirli  tligln 
data  have  been  developed.  These  proredures  have  been  applied  to  a  nmtiber 
oflateral  manoeuvres  for  wing  sweeps  of  72..')  to  16  degrees.  The  comparison 
of  model  predictions  with  flight  measurements  oflateral  response  reveal  the 
following  characteristics. 

•  A  very  good  degree  of  matching  between  the  model  calculated  values 
of  p.  r  and  i3  and  the  measured  values  is  possible. 

•  For  some  cases  the  yaw  rate  r  is  poorly  estimated.  This  is  believed  to 
be  due  to  inaccurate  C„  derivatives. 

•  When  the  yaw  rate  is  inaccurate,  the  accuracy  of  the  calculated  sideslip 
angle  J  is  significantly  reduced,  whilst  the  accuracy  of  the  roll  rate  p 
is  affected  only  to  a  minor  degree. 

•  The  derivative  is  often  over-estimated,  which  degrades  the  accu¬ 
racy  of  the  calculated  roll  response. 

•  Strong  coupling  exists  in  the  sideforce  equation  arising  from  resolving 
roll  inertia  terms  into  the  airpath  axes.  Errors  in  the  calculation  of 
angle  of  attack  can  lead  to  incorrect  resolution  of  these  inertia  terms. 

Results  of  a  preliminary  investigation  of  a  longitudinal  manoeuvre  have 
been  presented.  These  results  show  that  the  wind-tunnel  data  for  the  lon¬ 
gitudinal  aerodynamic  parameters  is  much  more  accurate  than  the  lateral 
data.  Hence  smaller  adjustments  to  the  data  will  be  needed. 

Further  investigation  into  the  calculation  of  the  yaw  rate  to  determine 
the  suspected  inaccurate  C„  derivatives,  is  recommended. 
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p,r,j3  calculated  seperately 


p,r,0  calculated  simultaneously 


Fig. 


- measured  data 

—  —  calculated  data 

1  Comparison  of  measured  and  calculated  data  of  a  lateral 
manoeuvre,  A  =  72.5°,  using  Flight  Trials  derivatives. 


p,  r,  /?  calculated  seperately 


p,T,P  calculated  simultaneously 


- measured  data 

■■■  calculated  data 

Fig.  3  Comparison  of  measured  and  calculated  data  of  a  lateral 
manoeuvre,  A  =  35®,  using  Flight  TrieJs  derivatives. 


r 


Fig.  6  Improvement  in  calculation  of  r  (seperately  from  other  variables) 
due  to  altering  C„  derivatives. 


q,a,9  calculated  seperately  9,o,^  calculated  simultaneously 

- measured  data 

- calculated  data 

Fig.  7  Comparison  of  measured  and  calculated  data  of  a  lonitudinal 
manoeuvre,  A  =  16°,  using  Database  derivatives. 


Fig.  10  Inaccurate  estimation  of  yaw  rat 
(Results  for  A  =  72.5,35.4:26”  ). 


Fig.  0  Example  of  close  matching  of 
measured  and  calculated  p,r,^ 
variables  calculated  simultaneously. 
(Results  for  A  =  50°  ) 


a 


Fig.  lib  Yaw  Tale  calculated  seperately 
with  altered. 


Fig.  12  Yaw  rate  calculated  fimultaneoualy  with 
Cng,Cn,,  altered  as  in  Fig.  11b 

■  calculated  data 
- measured  data 

Results  for  A  =  35°  lateral  manoeuvre 


Fig. 


13a  0,p  calculated  simultaneously 
with  correct  yaw  rate  as  input. 


~  ••  “•  =•“ 

Fig.  ISb  calculated  simultaneously 
with  calculated  yaw  rate  (shown 
in  Fig.  12)  as  input. 


d.  ud 


10  Cl 


time 

Fig.  14  Calculated  and  measured  with 

Cy  derivatives  equal  to  flight  trial  values. 


Fig.  16  FYS  term  with  Cy  derivatives 
equal  to  flight  trial  values. 


derivatives  equal  to  zero. 


Fig.  15b  Comparison  of  and  Fig.  17b  Calculated  and  measured  13  with 

TFYW  terms.  Cy  derivatives  zero.  (Compare  to  Fig.  14) 


Fi«.  20  Calculated  roll  rate  with  overestimated  value, 
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Fig.  21  Calculated  roll  rate  with  reduced  C/,,  value. 


■  ■  calculated  data 
- -  measured  data 
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